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Abstract 

Bose-Einstein correlations of both neutral and like-sign charged pion pairs are 
measured in a sample of 2 million hadronic Z decays collected with the L3 detector at 
LEP. The analysis is performed in the four-momentum difference range 300 MeV < 
Q <2 GeV. The radius of the neutral pion source is found to be smaller than that of 
charged pions. This result is in qualitative agreement with the string fragmentation 
model. 
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Introduction 



Bose- Einstein correlations (BEC), between identical bosons, have been extensively studied in 
hadronic final states produced in e+e~, ep, hadron-liadron and heavy- ion interactions The 
bosons studied are mainly charged pions @J^. Only rarely have neutral pions been studied 
and never before in e"'"e~ interactions. In this letter, we report a study of BEC of 7r° pairs in 
hadronic decays of the Z boson at LEP, and compare them to BEC of pairs of identically 
charged pions. 

BEC manifest themselves as an enhanced production of pairs of identical bosons which are 
close to one another in phase space. This can be studied in terms of the two-particle correlation 
function R2 @]: 

^2 Pl,P2 = 7 ^, 1) 

Po{Pl,P2) 

where Pbe{pi,P2) is the two-particle number density for identical bosons with four-momenta joi 
and P2, subject to Bose-Einstein symmetry. The reference distribution, Po{pi,P2), is the same 
density in the absence of Bose-Einstein symmetry. 

Assuming a static spherical boson source with a Gaussian density and a plane wave descrip- 
tion of the bosons, -^2(^15^2) is written as 

R2{Q) = Afil + aQ){l + Ae-«'^'), (2) 

where = —{pi— P2Y is the square of the four-momentum difference. The parameter R 
can be interpreted as the size of the boson source in the centre-of-mass system of the boson 
pair and a measurement of the correlation function R2 gives access to the source size. The 
parameter A is introduced to describe the fraction of effectively interfering pion pairs. In this 
analysis the normalization factor M{l + aQ) is added. It takes into account possible long-range 
momentum correlations, as well as possible differences in pion multiplicity in the data and 
reference samples, over the four-momentum difference range studied. 

The spherical shape of the boson source assumed here is a simplified picture. High statistics 
charged pion data at LEP revealed the source to be elongated PJTO[| . The present measurement 
has, however, no sensitivity to the shape of the source of neutral pions, due to its limited 
statistics. 

Several theoretical predictions exist for differences in BEC for pairs of bosons in the pion 
isospin triplet (Tr"*", tt", 7r°). From the string model |11] a smaller spatial emission region, i.e. 



a wider momentum correlation distribution is expected for tt'^tt^, than for vr^vr^. This follows 
from the break-up of the string into qq pairs, which forbids two equally charged pions to lie 
next to each other on the string, whereas two neutral pions can. The same effect is found when 
the probabilistic string break-up rule is interpreted as the square of a quantum mechanical 



amplitude From a quantum statistical approach to Bose-Einstein symmetry [jl4|, a 

small difference between vr^vr^ and tt^tt" correlation is expected. The size and shape of this 
difference is predicted to be similar to an expected Bose-Einstein correlation of Tr^Tr"^ pairs. It 
is theoretically uncertain whether BEC between unlike sign pions are also to be expected on 



the ground of isospin invariance |]15 



The main purpose of this letter is to measure the difference in size of the emission region of 
neutral and charged pions. In order to minimize systematic uncertainties on this difference, the 
procedures followed for the charged pions are kept as close as possible to those for the neutral 
pions. 
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Data and Monte Carlo Samples 



For this analysis, data collected with the L3 detector [0 during the 1994 and 1995 LEP runs 
are used. The analysis is based mainly on measurements from the high resolution electro- 
magnetic calorimeter and from the central tracking device. The data sample corresponds to 
an integrated luminosity of about 78 pb~ at centre-of-mass energies around ^/s = 91.2 GeV. 
From this sample, about 2 million hadronic events are selected, using energy deposits in the 
electromagnetic and hadronic calorimeters |T^. 

The Jetset generator [|1^] is used to study the detector response to hadronic events. Param- 
eters of the generator are tuned to give a good description of event and jet shapes of hadronic 
events measured in L3. The effects of Bose-Einstein symmetry are simulated with the LUBOEI 
routine [19]. The routine has two parameters, which have been chosen to obtain a reasonable 
description of L3 data. This ad-hoc model shifts boson momenta after the hadronization phase 
in such a way that the correlation function i?2(Q) for identical bosons is proportional to a 
constant plus a Gaussian as in Equation 0. The generated events are passed through a full 
detector simulation and are reconstructed and subjected to the same analysis procedure as 
the data. This Monte Carlo sample (Jetset-BE) contains about 7 million events. A control 
sample (Jetset-noBE) is also generated with Jetset but it has LUBOEI switched off and 
the generator parameters retuned. Significant differences are found in the tuned parameters in 
the two cases. The number of events in the control sample is approximately 2 million. Unless 
stated otherwise, the Jetset-BE sample is used throughout this letter. 



Neutral Pion Selection 

Neutral pions in hadronic events are reconstructed from photon pairs. Photon candidates are 
identified in the electromagnetic calorimeter as a cluster of at least two adjacent crystals. The 
clusters are required to be located in the central region of the detector, | cos(^ciuster)| < 0.73, 
and to be in the energy range 100 MeV < -^cluster < 6 GeV. Above 6 GeV, the two photons 
from a 7r° decay can no longer be distinguished as two separate clusters. 

Discrimination of clusters originating from photons or electrons from those due to other 
particles is based on the distribution of the energy over the crystals of the cluster. A good 



photon discrimination is achieved with a neural network based on this energy distribution ||2T 
To reject clusters due to charged particles, a minimum distance between the cluster and the 
extrapolation of the closest track of 30 mm is required. This corresponds to 1.5 times the size 
of the front face of a crystal. 

Pairs of photon candidates within an event are used to reconstruct vr^'s. The distribution of 
invariant mass of photon pairs shows a peak around the tt" mass, above a smooth background. 
These two components are extracted by a fit to the mass spectrum. The background is described 
by a Chebyshev polynomial of third order, /bg. The 7r° peak is parameterized by a Gaussian 
function with exponential tails, which is continuous, and smooth in the first derivative: 



exp I — 1 exp I j if m.y^ — < —aa 



'111 



exp \ -- ^^^2 j if /5cr > ''^77 — m-n > —aa . (3) 

exp (y2^ exp ^_ /^(™77^-™-vr) j •£ _ ^ 

Here, m^^ is the two-photon invariant mass, indicates the peak position, and a is the width 
of the Gaussian peak. The parameters a and (3 determine the values of m^^ where the Gaussian 
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changes into an exponential. The exponential tail at low mass is needed to describe the presence 
of converted photons in cluster pairs. If a photon converts, e.g., in the outer wall of the tracker, 
it can still be selected as a photon, although some of the original photon energy is lost. The 
high mass tail accounts for an overestimate of the cluster energy due to another cluster nearby. 

An example of the fit result is given in Figure 0. The sum of the two functions /bg and /tt 
describes the distribution well. Also, the shape and the size of the 7r° peak in the data and 
Monte Carlo agree. The mass resolution as determined from the fit is about 7.4 MeV. 

The photon pair is then selected as a 7r° candidate if it has an energy 200 MeV < E^^ < 
6 GeV. For a mass window 120 MeV < m^^ < 150 MeV, a total of 1.3 million vr^'s is selected 
in data. The 7r° purity of the candidate sample is of the order of 54% and the 7r° selection 
efficiency is approximately 17%. About half of the background is combinatorial, i.e., photons 
from different tt^'s. The other half consists of pairs where one or both of the photons do not 
come from a tt" decay. 

Charged Pion Selection 

Charged pions are detected as tracks in the central tracker. They are selected in the same 
kinematic range as neutral pions: | cos(6'track)| < 0.73 and 200 MeV < Etmck < 6 GeV, where 
the energy is calculated from the track momentum, assuming the tt^ mass. In addition, at least 
35 out of 62 possible wire hits are required in the track fit, and the number of wires between 
the first and the last must be at least 50. Furthermore, the track must have at least one hit 
in the inner part of the tracker, and the distance of closest approach to the e"'"e~ vertex in 
the plane transverse to the beam is required to be less than 5 mm. Finally a high resolution 
measurement of the polar angle is demanded. Charged pions are analyzed in the 1995 data 
only, in which 4.1 million tracks are selected. 



Pion Pair Analysis 
Neutral Pions 

After the neutral pion selection, 7r° candidates within an event are paired, requiring that no 
cluster is common to the two candidates, and their four-momentum difference Q is calculated. 
The 7r°7r° component of the Q distribution is estimated by a fit to the two-dimensional mass 
spectrum for every bin in Q. An example of these mass distributions is shown in Figure for 
the bin 0.48 < Q < 0.52 GeV. The various contributions are clearly visible: non-yr'' pairs give 
the smooth background, 7r° with non-vr^ pairs give the two "ridges" in the peak regions, 7r° 
pairs give part of the peak in the centre of the plot, the other part being caused by the sum of 
the 7r° with non-7r° ridges. 

This two-dimensional distribution is derived from the product of two one-dimensional mass 
distributions PT[ |: 



/2d(mi,m2) = A^^f^{mi)f^{m2) 

+ ^7rbg[/7r(W-l)/bg(m2) + /bg("T-l)/7r(?^2)] (4) 
+ ^bgbg/bgbg("li, 777,2), 

where the first term describes the part, the second term describes the 7r° with non-7r° pair 
ridges and the third term is the non-7r° pair background. The number of pairs in the mass 
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window 120 MeV < mi 2 < 150 MeV is related to the parameter A^^. The functions f^, and 
/bg have the same functional form as described before, /bgbg follows from the product of two 
Chebyshev polynomials of third order and is required to be symmetric in the two masses. All 
18 parameters of /2d are left free in the fit. 

The result of a binned maximum likelihood fit for the mass distribution in Figure ^ is shown 
in Figures ^o-d. In this representative example, the is 4737 for 4606 degrees of freedom, 
which corresponds to a 9% confidence level. 

Figure Qa presents the Q distribution for pairs in the mass window 120 MeV < mi^2 < 
150 MeV, obtained from the values of the parameter A.,^^^ from two-dimensional mass fits of 
Equation ^to both data and Monte Carlo. Some deviations between data and Monte Carlo are 
caused by the imperfect modeling of BEC. The efficiency to select a 7r° pair in an event ranges 
from about 1% at Q = 300 MeV to 4% at Q = 2 GeV. 

Bins in Q below 300 MeV are not used for the rest of the analysis for two reasons. First, 
the efficiency estimate depends strongly on the BEC modeling in the generator, in the region of 
small Q. This occurs because the BEC modeling moves identical pions closer together, which 
lowers the detection efficiency. Secondly, the four- momentum difference of any pair of tt^'s from 
rj —>■ TT^Tr^Tr" decays is kinematically constrained to have Q < 311.7 MeV, and in that Q-range, 
more than 20% of all vr" pairs originate from this decay. A rejection of the small Q region thus 
avoids systematic uncertainties due to the simulation of the rj multiplicity. 

Charged Pions 

The distribution of four-momentum difference of equally charged pion pairs is obtained by 
calculating Q for pairs of tracks selected within an event and with the same charge. This raw 
TT^TT^ spectrum is corrected bin-by-bin for both pion purity and efficiency using the Monte 
Carlo simulation. The uncorrected distribution is shown in Figure ^3. Compared to the n^n^ 
case, smaller deviations are observed between the raw spectrum in data and Monte Carlo. 
These deviations are due to the imperfect modeling of BEC in the Monte Carlo. 

Results 
Neutral Pions 

To obtain the final correlation function R2{Q), the Q distribution of vr^ pairs in the data. 
Figure is corrected for selection efficiencies. The efficiency is defined as the number of 
selected 7r° pairs in Monte Carlo events (Jetset-BE) divided by the number of generated 7r° 
pairs in the same events, where the generated pions have to be in the same kinematic range as 
the selected pions. This definition includes an acceptance correction for those tt^'s which cannot 
be selected kinematically. In this way, the correlations of pairs can directly be compared 
to those of charged pion pairs. The reference distribution po{Q) is calculated from a Jetset- 
NoBE sample at generator level, where pions are taken in the same kinematic range as in 
the definition of the selection efficiency. We choose this reference distribution rather than the 
distribution for tt^tt^ because of the uncertainty concerning BEC between unlike sign pions 
mentioned in the introduction. The correlation distribution R2{Q) is then the ratio of the 
corrected data spectrum to the reference spectrum. 

The distribution of R2{Q) is displayed in Figure An enhancement at low Q values, 
expected from Bose- Einstein symmetry, is clearly visible. The function R2{Q) from Equation ^ 
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is fitted to this ratio in the interval 300 MeV < Q < 2 GeV. Extending the fit to lower values of 
Q results in consistent values of the parameters but with much larger systematic uncertainties. 
The overall normalization J\f is determined from the integrals of i?2(Q) and the fit function; 
the only free parameters are A, R and a. In this fit, the is 46.1 for 40 degrees of freedom, 
corresponding to a 23% confidence level. 

The systematic uncertainty on the result due to the 7r° selection is determined by varying 
the photon selection cuts and by changing the size of the tc^ mass window by ±10 MeV. The 
^0^0 j^g^gg f[f^ Qf Equation ^ is tested by varying the fit range by ±12.5 MeV. The uncertainty 
due to the modeling of Bose-Einstein correlations in the Monte Carlo generator is taken into 
account by using the control sample Jetset-noBE in the efficiency correction procedure. In 
addition, the influence on the final result of the agreement between data and Monte Carlo 
of distributions relevant to the photon and 7r° selection, such as neural network output and 
energy and polar angle of 7r°'s, is studied. Finally, the binning in Q is varied. The systematic 
uncertainty on the result due to each of the sources, is assigned as half the maximum deviation. 
A summary is given in Table |I|. The total systematic uncertainty is calculated as the quadratic 
sum of these uncertainties. 

Charged Pions 

The final correlation function i?2(Q) for tt^tt^ is calculated in a similar way as that for n^n^. 
The TT^TT^ data distribution. Figure §3, is corrected for purity and efficiency. As for vr'^vr'^, the 
efficiency is calculated for generated pions in the same kinematic range as the selected pions. 
The reference distribution po{Q) and the correlation distribution i?2(Q) are obtained in the 
same way as for tt^tt". 

The correlation function for tt'^tt^ is shown in Figure Qd. Due to the higher selection 
efficiency for charged pions as compared to neutral pions, the significance of the low Q value 
enhancement is much larger. As for tt^tt^, the function defined in Equation ^ is fitted to the 
final distribution. In this fit, the is 42.6 for 40 degrees of freedom, corresponding to a 36% 
confidence level. 

The systematic uncertainty on the result due to the track selection is determined by varying 
the requirements on number of hits, distance of closest approach and polar angle determination. 
As in the n^n'^ case, the uncertainty on the modeling of Bose-Einstein correlations in the Monte 
Carlo generator is obtained by using the control sample Jetset-noBE in the analysis. Finally, 
the binning in Q is varied. The systematic uncertainties are attributed as in the ir^n^ case, 
and are summarized in Table 0. The total systematic uncertainty is calculated as the quadratic 
sum of these uncertainties. 

Comparison 

The final values for the strength of the correlation A and the corresponding radii of the boson 
sources R are given in Table 

Due to the lower efficiency of the n^Tc^ selection, the statistical uncertainty on the tt^tt^ 
result is larger than the statistical uncertainty on the tt^tt^ result. Within these uncertainties, 
the data indicate both a weaker correlation and a smaller source radius for tt^tt^. The weakness 
of the 7r°7r° correlation can be partly explained by the bigger contribution of resonance decays 
to the Q spectrum. The difference of the source sizes is 

- Roo = 0.150 ± 0.075 (stat.) ± 0.068 (syst.) fm, (5) 
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where R±± and i?oo indicate the value of R for tt^'^tt^ and tt^tt*^, respectively. In this difference, 
the systematic uncertainties due to the modeling of Bose-Einstein correlations and the binning 
in Q are taken to be correlated between the two samples. The smaller radius found for 7r°7r° is 
in qualitative agreement with the predictions of the string model. 
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Source 


AA 


AR (fm) 


Aa 


Photon selection 


0.020 


0.052 


0.009 


Mass window 


0.011 


0.008 


0.007 


2D fit range 


0.056 


0.019 


0.036 


MC modeling 


0.037 


0.035 


0.020 


Data-MC agreement 


0.012 


0.018 


0.034 


Q-binning 


0.004 


0.013 


0.003 


Total 


0.072 


0.070 


0.055 



Table 1: Systematic uncertainties on A, R and a for the 7r°7r° data sample. 



Source 


AA AR (fm) Aa 


Track selection 
MC modeling 
Q-binning 


0.011 0.009 0.009 
0.022 0.003 0.002 
0.001 0.001 0.001 


Total 


0.025 0.010 0.009 



Table 2: Systematic uncertainties on A, R and a for the n^n^ data sample. 



Sample 


A R (fm) a 




0.155 ±0.054 ±0.072 0.309 ± 0.074 ± 0.070 0.021 ± 0.034 ± 0.055 
0.286 ± 0.008 ± 0.025 0.459 ± 0.010 ± 0.010 0.015 ± 0.003 ± 0.009 



Table 3: Values for A, R and a, for both the n^ir^ and the tt tt data samples. The first 
uncertainty is statistical, the second systematic. 
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Figure 1: Distribution of (a) the two-photon invariant mass m^^ for data together with the 
fit result and (b) the 7r° signal as obtained from fits to data and to Monte Carlo. The arrows 
indicate the mass selection window. 
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Figure 2: (a) Two-dimensional distribution of the mass of tt^ pair candidates with a four- 
momentum difference in the ransr e 0.48 <Q < 0.52 GeV, in bins of 2.5 x 2.5 MeV^ Result of 
the fit for (b) tt*^ pairs, (c) 7r° with non-7r° pairs, (d) non-Tr" pairs. 
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Figure 3: Data and Monte Carlo distribution of the four- momentum difference of (a) pairs of 
TT^'s, as obtained from fits of Equation ^, and (b) pairs of tt^ candidates. 
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Figure 4: Distribution of i?2(Q) for (a) ti^ti^ and (b) tt'^tt^, and results of the fits. The 
points indicate the data, the full line corresponds to the fit result and the dashed line is the 
normalization factor M{1 + oiQ). 
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